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Abstract 



The phase-ordering kinetics of a bulk uniaxial nematic liquid crystal is ad- 
dressed using techniques that have been successfully applied to describe or- 
dering in the 0{n) model. The method involves constructing an appropri- 
ate mapping between the order-parameter tensor and a Gaussian auxiliary 
field. The mapping accounts both for the geometry of the director about the 
dominant charge 1/2 string defects and biaxiality near the string cores. At 
late-times t following a quench, there exists a scaling regime where the bulk 
nematic liquid crystal and the three-dimensional 0(2) model are found to be 
isomorphic, within the Gaussian approximation. As a consequence, the scal- 
ing function for order-parameter correlations in the nematic liquid crystal is 
exactly that of the 0(2) model, and the length characteristic of the strings 
grows as t^^^. These results are in accord with experiment and simulation. 
Related models dealing with thin films and monopole defects in the bulk are 
presented and discussed. 

PACS numbers: 05.70.Ln, 61.30.-v, 64.60.Cn 



1 



Typeset using REVTgX 



I. INTRODUCTION 



Most phase-ordering systems studied to date support only one type of topologically 
stable defect species One example is the 0{n) model with an n-component vector 

order-parameter. In three spatial dimensions, the defects formed at the quench are line-like 
strings for n = 2, and point-like monopoles for n = 3. Phase-ordering in bulk uniaxial 
nematic liquid crystals (nematics) provides the simplest scenario in which two defect species 
- monopoles and strings - are topologically stable. The stability of monopoles derives from 
the 0(3) symmetry of the nematic director n{f, t). The additional invariance under the local 
inversion n(r, t) —n{f, t) allows the nematic to support stable charge 1/2 disclination lines 
(strings) 0. The issue of which defect species dominates the dynamics in bulk nematics 
at late times t following a quench has recently been settled. Cell-dynamical simulations 
using spin models of bulk nematics have computed the order-parameter correlation 
function and found it to be indistinguishable from that of the 0(2) model, and consistent 
with a string-dominated late-time scaling regime. Experiments by Chuang et al. |]^ directly 
imaged the bulk nematic, revealing an intricate, evolving defect tangle. While both types of 
defect were observed, the strings dominated at late times. The length scale Ls characterizing 
the typical separation of the strings was seen to grow as Lg ~ t^^^ while the average line 
density of string {if) decayed like {rf) ~ ~ t^^ . The study of ordering in nematics is also 
of interest to cosmologists since similar processes involving cosmic string and monopole 
evolution, thought to occur in the early universe, may be responsible for structure formation. 

In this paper a theory is presented that describes the dominant scaling behaviour of the 
bulk nematic in terms of a string-dominated late-time regime. Generalizing a successful 
method used to treat the ordering kinetics of the 0(n) model, the nematic order-parameter 
tensor is mapped onto a two-component Gaussian auxiliary field 0. The string defects 
explicitly appear in the construction of the mapping. As discussed below, this approach has 
several advantages over an earlier, semi-numerical theory by Bray et al. [|l^ . 



The auxiliary field approach is first applied to the straightforward case of phase-ordering 
in nematic films containing charge 1/2 vortices, which have been studied in simulations 
and experiments |TT]]. As in the bulk nematic, the mapping is constructed to account for the 
rotation of the director by only vr about the core of the defect. Once this is done the theory 
reveals that phase-ordering in the nematic film is equivalent to phase ordering in the two- 
dimensional 0(2) model examined previously 0. This is not surprising since the two systems 
are known to be isomorphic |]5|JT^. Constructing a theory for the bulk nematic is more 
challenging since the order-parameter tensor must include a biaxial piece near the core of the 
string. In the earlier theory of Bray et al. |jlO| this point was not addressed since there they 
used a "hard-spin" approximation for the dynamics of the nematic. However, the necessity 
of a having a biaxial core region when treating the full equations has been noted in the 
numerical work of Schopohl and Sluckin |]TB[ on bulk nematic string defects in equilibrium. 
The present theory successfully incorporates biaxiality and clarifies the role that it plays in 
the coarsening of the bulk nematic. The theory recovers the growing length Lg ~ t^^^ seen in 
simulations and experiments 0]. In the scaling regime, the order-parameter correlation 
function for the bulk nematic is found to be exact/?/ that of the three-dimensional 0(2) model 
0, in excellent agreement with simulations [^] (Fig. [1|). Although the theoretical results of 



Bray et al. |10] suggested agreement between the correlation function for the bulk nematic 
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and the 0(2) model they were unable to demonstrate an exact equivalence since their theory 
was not based on a mapping that explicitly contained strings. The major accomplishment 
of this work is to analytically demonstrate the isomorphism between the dynamics of the 
bulk nematic and the dynamics of the three-dimensional 0(2) model, within the Gaussian 
approximation. Through this isomorphism, the well-developed theory for the 0(2) model 



T^JT5|| can be applied directly to the nematic. In particular, this theory predicts that 



the average line density of string decays as {rj) ~ ~ t [0,^, in accord with the 
experiments of Chuang et al. 0. 

Although strings are generically present in bulk nematics, certain choices of experimental 
setup and sample material will produce copious amounts of monopoles at the quench ^\ . 
The theory of Bray et al. [|l^ is unable to address these experiments since in that theory 
there is no signature for monopoles. However, within the framework presented below it is 
relatively straightforward to develop a theory of nematics in which monopoles appear. In this 
theory the order-parameter correlation function is found to be similar to (but not exactly) 
that for the three-dimensional 0(3) model [0] (Fig. 0). The characteristic monopole spacing 
Lm grows as Lm ~ ^ ^ and leads to a decaying average monopole density (n) 
Experiments that examine monopole-antimonopole annihilation in isolation from strings 
suggest that these growth laws should hold. However, experiments [0 also reveal that the 
average monopole density decays more rapidly in the presence of strings, with (n) ~ t~^. It 
appears that in order to account for this observation the theory presented here should be 
extended to consider the interactions between strings and monopoles ||16| . 



II. MODELS 

In this section the 0(n) model and the Landau-de Gennes model of nematics are dis- 
cussed. Since the former model is used as a guide in the treatment of the latter, the theory 
for ordering kinetics in the 0(n) model is also reviewed. Initially, the structural features 
common to both models are emphasized. In later sections, the technical details specific to 
the ordering of nematics will be discussed. 



A. The 0(n) model 

In the 0(n) model the evolution of the non-conserved, n-component order-parameter 
field ifj is governed by the time-dependent Ginzburg-Landau (TDGL) equation 

dt 5^ ' ^ ' ' 

The free energy F[il)\ has the form 

F[^]= j d'r[]^{y^Y + Vm (2.2) 

where the potential V^(V'); expressed in terms oiip = is 0{n) symmetric with a degener- 
ate ground state at non-zero = ipQ. In this model, as with the nematic liquid crystal, the 
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disordered high temperature initial state is rendered unstable by a quench to a low tempera- 
ture where the usual noise term on the right-hand side of (|2.1| ) can be ignored. Substitution 
of ( p.2| ) into produces the explicit equation of motion 



^.VV-^. (2.3) 
dt 



The evolution induced by (|2.3|) causes ip to order and assume a distribution that is far 
from Gaussian. To make analytic progress it is by now standard to introduce a mapping 

if=a{m) (2.4) 

between the physical field and an ra-component auxiliary field fh with analytically tractable 
statistics. The mapping a is chosen to reflect the defect structure in the system and satisfies 
the Euler-Lagrange equation for a defect in equilibrium 



As shown below, ( |2.5| ) is also instrumental in treating the non-linear potential term in 
Defects correspond to the non-uniform solutions of ( |2.5| ) which match on to the uniform 
solution far from the defect core. Since only the lowest-energy defects, those with unit 
topological charge, will survive until late-times the relevant solutions to (|2.5| ) will be of the 
form P] 

(T(m) = A{m)'m (2.6) 

where m = |m| and m = m/m. Thus the magnitude of m represents the distance away from 
a defect core and its orientation corresponds to the orientation of the order parameter field at 
that point. This geometrical interpretation will later be exploited when the generalization of 
( p.5|) is used to choose the appropriate mapping, analogous to ( p.6|) , for string defects in the 



nematic liquid crystal. The magnitude of m grows as the characteristic defect separation, 



L{t), becoming large in the late-time, scaling regime. Inserting (|2.6|) into (|2.5| ) gives an 
equation for A, the order-parameter profile around a defect 

- - ^{A) = 0. (2.7) 

oA 

For small m an analysis of ( |2.7| ) yields the linear dependence A(m) ~ m, characteristic of a 
unit charge defect [|18|. For large m the amplitude A approaches its ordered value A = ipQ 
algebraically, which is a feature common to both the 0{n) model and the nematic. 
The order-parameter correlation function is 

C(f,t) = {a{r,t)-a{Q,t))=^l{m{r,t)-m{(),t)) (2.8) 

where the last equality holds for late-times and to leading order in 1/L. To evaluate the 
last average in (|2.8|) we choose m to be a Gaussian field with zero mean. This Gaussian 
approximation forms the basis of almost all present analytical treatments of phase-ordering 
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problems, and has had much quantitative success in describing the correlations in these 
systems Theories where m is a non-Gaussian field also exist |19;20]. In the Gaussian 

approximation the order-parameter correlation function ( |2.8| ) can be related to the normal- 
ized auxiliary field correlation function /, defined as 



f{r,t) 



{m{f^ t) ■ m(0, t)) 



The relation is P,|21 



([m(0,tP) 
C{r,t)=i^l:F{r,t) 



with 



2-K 



I n+l 
2' 2 



1 1 n 



2' 2' 2 



(2.9) 



(2.10) 



(2.11) 



where B is the beta function and F is the hypergeometric function. In the late-time scaling 
regime the functions T and / can be expressed solely in terms of the scaled length x = r/L{t) 
so that JF = jF(a;). In this regime the equation of motion (|2.3| ) can be written as non-linear 
scaling equation for JF 



TT 



d 



(2.12) 



In the derivation of ( |2.12| ) the relation ( p.5|) is used to replace the potential term in ( p^ 
and then the Gaussian identity 



([V^a(m(f,t))]-a(m(0,t))) 



nf{f,t) 



d 



([m(0,t)]2) dfif,t) 



{a{m{f,t)) ■ a{m{0,t))) (2.13) 



is used to get the last term on the left-hand side of ( p.l2| ). The constant fi enters through 
the definition of the scaling length L: 



L\t) ^ 



2nfi 



At. 



(2.14) 



This is the well-known p|,^ growth law L ~ t^/^ for phase-ordering in non-conserved vector 
systems. 

Since the auxiliary field fh is smooth [jl5[, / is analytic for small-x. This implies, through 
an examination of ( p. 12 ) in d spatial dimensions, that for small- a; J-' behaves like 



TT 



jTfx) = 1 + — -x^ In x + C(x^l 



for n = 2 and 



^(x) 



TT 



TC 



2nd 3jj{d+l)V2nd 



X 



0{x^) 



(2.15) 



(2.16) 



for n = 3, the cases relevant to this paper. The non-analytic terms in refiect the short- 
distance singularities in the order parameter field produced by the defects, and lead to the 
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Porod's law power law decay of the structure factor at large wavenumber. The In x 
term in ( p.l5|) is characteristic of string (or vortex) defects while the term in ( 2.16|) is 
due to monopole defects. For large x both and / decay rapidly to zero. The eigenvalue 
/i is determined numerically by matching the short- and long-distance behaviours of the 
solution of ( p.l2| ). In this way the auxiliary field correlation function / is determined self- 
consistently along with JF. In contrast, there is no such self-consistency in theories based 
on the Ohta-Jasnow-Kawasaki (OJK) approximation |T^J2^. Values of ^ at various n and 
d for the 0{n) model have been determined 0. The scaling functions JF of this theory are 
in excellent agreement with the results of simulations 



B. Nematic Liquid Crystals 

The order-parameter for a bulk nematic liquid crystal is a traceless, symmetric, 3x3 
tensor Qq,/3, which measures the anisotropy of physical observables in the nematic phase. 
The tensor Qap has the general form 

1 1 

Qap = A[hahp - -5„/3] + ^B[gc,gp - hahp]. (2.17) 

The unit 3-vectors h, g and h form an orthonormal triad. The amplitudes A and B are 
chosen to be non-negative. A is a measure of the degree of uniaxial order in the liquid 
crystal; it is zero in the isotropic phase and non-zero in the nematic phase. Biaxiality in the 
liquid crystal is measured by i?. In the uniaxial nematic phase B is zero everywhere except 
near the string cores. The description of nematics in terms of Qai3 reduces to the Frank 
continuum theory of elasticity in terms of a director n{f, t) when A is set to its ordered 
value and 5 = 0. In the phase-ordering scenario, where defects occur, all of A, B, fi, g and 
h are space and time dependent. 

In the tensor formulation the director, which measures the average local molecular ori- 
entation in the nematic, is the unit eigenvector of Qais which corresponds to the largest 
eigenvalue. The unit eigenvectors and associated eigenvalues of Q^fs are: 

h ^ |y4 

g ^ -liA-B) (2.18) 
h ^ -1{A + B) 

Since the nematic is uniaxial, B < 3A and the director can be identified with n. The tensor 
formulation respects the full RP^ symmetry of the uniaxial nematic since physical quantities, 
such as correlations, are written in terms Qa/s which is invariant under the local inversion 
h{r,t) — >■ —h{r,t). At a string core B = 3 A > and the eigensubspace corresponding to 
the largest eigenvalue 2 A/3 is two- fold degenerate. Thus in the plane perpendicular to h, 
the tangent to the string, the orientation of the director is ambiguous. At the isotropic core 
of a monopole A = B = and all three eigenvalues of Q^/? are degenerate so the orientation 
of the director is completely unspecified. 

The dynamics of the nematic is governed by the TDGL equation for Qa/B 

dtQap = + KpTrQ (2.19) 
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with the Lagrange multipher Xap included to enforce the traceless condition. The Landau-de 
Gennes free energy is 

F[Q]= Jd\[^{VQ)' + Vm (2.20) 

with the potential 

V{Q) = -^TrQ' - ^TrQ' + \iTrQY. (2.21) 

The coefficient of the quadratic term in ( p.21|) is chosen to be negative so that the bulk 
isotropic phase is unstable towards nematic ordering. The gradient term in (p.20|) is written 
within the equal-constant approximation |2^. Substitution of the form ( |2.17| ) in ( p.21| ) 
results in a useful expression for the potential as a function of A and B: 

V{A, B) = --A^ - —A" + -A^ - —B^ + —B^ + —[AB^ + A^B\ (2.22) 
^ ' ^ 9 27 9 27 81 27^ ^ ^ ' 

A contour plot of V{A, B) for A > 0, i? > is shown in Fig. ^ There is a global 
isotropic maximum at {A,B) = (0,0), a uniaxial minimum at {A,B) = (1,0), and a saddle 
at {A,B) = (1/4,3/4). The minimum represents the bulk nematic phase, the isotropic 
maximum corresponds to the monopole core, and the saddle, with B = 3A, signifies the 
string core. 

Substituting ( |2.20| ) into ( p.l9| ) and using TrQ = to calculate the Lagrange multiplier 
gives an explicit form for the equation of motion 

dtQa/3 = — Pajd (2.23) 

with the non-linear piece given by 

Pap = {TrQ' - l)Qap - [Q'U + ^-fTrQ\ (2.24) 
Static solutions to ( [^.23[ ) satisfy the Euler-Lagrange equation 

y^QaP = Pats. (2.25) 

The order-parameter correlation function is defined as 

C(r, t) = Af{Tr[6Q{f, t)6Q{0, t)]) (2.26) 

where Af = {{Tr[6Q{0,t)6Q{0,t)]))^^ is a normalization factor and 6Qai3 = Qap — {Qap)- 
Both Af and {Qap) are constant at leading order in 1/t. Thus, at late-times, equation ( |2.23| ) 
can be written as an equation for the evolution of order-parameter correlations 

UtC{r, t) = V2C(f, t) - Af{Tr[P{r, t)6Q{0, t)]). (2.27) 

Later, through a development that closely parallels that previously given for the 0{n) model, 
it will be shown how ( |2.25| ) and ( |2.27|) lead to a scaling equation for order-parameter corre- 
lations in the nematic. 
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III. STRING DEFECTS IN THE NEMATIC 



At late-times the dominant defects in the bulk nematic are strings with topological charge 
1/2. Many of the main features of phase-ordering in the bulk nematic are described by the 
model containing strings which is presented in Sec. III.B below. 



A. Vortices in thin films 

To begin, a model applicable to nematic thin films where the director is constrained to lie 
in a plane without breaking the h —>■ —n symmetry is examined. By restricting the director 
to a plane, the intricacies of how to map the order-parameter tensor onto an auxiliary field 
when the director rotates by only vr about the vortex can be demonstrated, without the 
additional complication of biaxiality that appears near the string core in bulk samples. 

For a uniaxial thin film nematic the order-parameter is a 2 x 2 traceless symmetric tensor 



(3.1) 



where h is the two-component director. In analogy to the theory of the 0(2) model, the 
defects are incorporated through a mapping of the order-parameter tensor onto a two- 
component auxiliary field. The only defect species present at late-times are charge 1/2 
point vortices with property that the director rotates by only vr around the vortex. This 
property is essential in constructing the mapping. Consider a charge 1/2 vortex at the origin 
with the typical director configuration 



h = cos -(p X + sin -0 y 



(3.2) 



where is the polar angle in the x — y plane. For future convenience we write the radial 
vector in the x — y plane as s and define angles in terms of s through 



[Sl,S2j 



COS0, sm( 



With the definitions ( p.2|) and (|3.3| ) the order-parameter tensor (|3.1|) is [26 

A{s) 



Q 



Sl S2 

h -Sl 



(3.3) 



(3.4) 



where s = This form for QajS, analogous to the mapping (2^) for the 0{n) model, is a 
solution to the Euler-Lagrange equation ( |2.25| ) written in terms of s 



VlQaP = Pa0 (3.5) 

where Pq/3 has a slightly modified definition from Pap ( p.24| ) because Qa/s is a 2 x 2 tensor: 



Pr 



a/3 



(Trg2 _ i/3)Q„^ - [Q\p + ^TrQ^ 



(3.6) 



Substituting ( p.4|) into (|3.5| ) results in an equation for the amplitude A: 
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From (|2.21| ) and (|3.4|) the potential U is given by 



U{A) = -^A' + ^A\ (3.8) 

An examination of (|3.7| ) at small s has A ~ s, indicative of charge 1/2 vortices ||TB[. At 
large s the amplitude A algebraically approaches its ordered value A = ^2/3. 

To treat many such vortices in a phase-ordering context, s in ( |3.4| ) is taken to be a 
Gaussian field s{f, t) with zero mean. As in the 0(2) model, s represents the distance to the 
nearest vortex, growing as the characteristic vortex spacing L^(t) at late-times. However, 
unlike in the 0(2) model, the director is not mapped directly onto s - a 27r rotation of s 
about a vortex corresponds to a rotation of the director by tt. 

At late-times, the amplitude A approaches it's ordered value and from the definition 
( p.26| ) and equation (|3.4|) order-parameter correlation function is seen to be 

C{f,t) = {s{f,t)-s{0,t)) (3.9) 

to leading order in L~^. This is just the 0(2) correlation function (|2.8|), and is related 
to /, the correlation function for the auxiliary field s defined in analogy to ( p.9|) , through 
( [^.l(j| ) and ( ^.11| ) for n = 2. In the scaling regime the equation of motion ( |2.27| ) for C{r,t) 
becomes equation (|2.12|) for the 0(2) scaling function J-", expressed in terms of the scaled 



length X = r/L^it). The length has the same definition as the length L in ( p.l4| ), with 
m replaced by s. The path from ( p.27| ) to ( |2.12|) is similar to that taken in the 0(2) case 
P]. The Euler-Lagrange equation (|3.5|) is used to replace occuring in the last term in 
( p.27| ). The resulting expression is evaluated using the Gaussian identity 



(Tr[V^Q(.-'(f,t))Q(.-'(0,t))]) = - ^[^(o^;^^ qj^^^ ^) (TrQ[(.-(r, t))Q(.-(0, t))]), (3.10) 



analogous to (|2.13|) , and produces the last term on the left-hand side of (|2.12|) . 



Thus the scaling function for the order-parameter correlations and the growth law 
L^{t) ~ t^/^ for the nematic thin film are exactly those of the two-dimensional 0(2) model. 
This correspondence, seen in simulations, can be simply understood as a consequence of 
the mapping of variables (p ^ 2(j) between the two models [@,|l3]. This isomorphism is 



relevant to experimental efforts that use constrained nematics to study coarsening in the 



two-dimensional 0(2) models |^ since it indicates that the existence of the local n —n 



symmetry does not affect the leading order dynamics in the scaling regime. 



B. Strings in the bulk nematic 

In addition to the complication of a having director configuration with a charge 1/2 
geometry, strings in a bulk nematic have a biaxial core. The form ( p.l7| ) for Qap contains 
the biaxiality that is required if an analytical solution to ( |2.25| ) is to be found. String defects 
enter the theory through the mapping of ( p.lTD onto a two-component auxiliary field. To 
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motivate the form for the mapping consider the geometry of the director field around a 
charge 1/2 string defect oriented along the z axis. Since locally the coordinate system can 
always be chosen so that the string has this geometry the following development is quite 
general. The director n is still given by (|3.2| ). The other members of the orthonormal triad 
in (|2T7|) are 



9 = - 

h = z 



sin -0 X + cos -0 y 



(3.11) 
(3.12) 



With the notation ( p.3|) for the radial vector s in the x — y plane, the order-parameter tensor 
( p.lTp becomes 



Q 



A{s) 











S2 
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1 - Si -S2 
-S2 1 + Si 

0-2 



(3.13) 



This form for Q^p is a solution of (|2.25|) written in terms of s, 



(3.14) 



provided that 



dA 

4 I jSA-B) dV 
3^^^ +3^^ ^dB 



0, 



(3.15) 
(3.16) 



where V{A, B) is given in ( p.22|) . Note that equations ( p.l5| ) and (|3.16|) would be inconsistent 
had a uniaxial ansatz {B = 0) been assumed at the outset. For the potential ( |2.22| ) these 
equations are degenerate [^] and reduce to a single equation for A after the identification 
B = 1-A: 



0. 



(3.17) 



At small s the solution to (|3.17|) is 



A = - + CS 
4 

B = --cs 
4 



(3.18) 
(3.19) 



where c is a constant, determined numerically. At large s the solution of ( p.l7| ) takes the 
form 



A 
B 



3 

4^ 



3 

4^2' 



(3.20) 
(3.21) 
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As expected, the mapping ( p.l3| ) connects the biaxial saddle point on the potential surface 
V{A,B), representing the string core, to the uniaxial nematic minimum away from the 
string (see Fig. |^). The linear behaviour in ( |3.18| ) and ( |3.19| ) at small s is that expected for 
charge 1/2 strings in the nematic. Both the linear behaviour near the core and the algebraic 
relaxation ( p.20| - |3l21| ) to the bulk uniaxial state are seen in the numerical results of [|1^]. 

Once again, to examine the statistical properties of the string defect tangle, s is a taken to 
be a Gaussian auxiliary field with zero mean. The magnitude s grows as the characteristic 
string separation Ls{t). Therefore, at late-times, s is large and the biaxial piece of Qai3, 
with an amplitude B given by ( |3.21 ), is suppressed. This is physically reasonable since 
biaxiality occurs on length scales around the core size while the late-time scaling properties 
are dominated by physics at the much larger scale of Ls{t). At late-times, when A ^ 1, the 
definition ( |2.26| ) and the mapping ( p.l3| ) show that the order-parameter correlation function 
reduces to 



C{f,t) = {sir,t) ■ siO,t)) (3.22) 

which is the 0(2) correlation function ( p.8|) . As before, C{f,t) is related to f{f^t), the 
normalized correlation function for the auxiliary field s, by relations ( ^.1(J| ) and ( |2.11| ) with 
n = 2. 

The dynamical equation ( p. 271) for C{f,t) reduces, in the scaling regime, to (|2.12|) for 
from the three-dimensional 0(2) model. Note that the spatial dimensionality enters 
through the Laplacian operator in ( p.l2| ). The scaled length in this case is a; = r/Ls{t) with 
Lg defined as L in ( |2.14| ). The derivation of this correspondence parallels the steps taken 
in the 0(2) model that lead to ( p.l2| ). The Euler-Lagrange equation ( p.l4| ) enables the 
non-hnear quantity Pap occuring in the last term of ( |2.27| ) to be replaced by V^Qa/?- The 
resulting average is then evaluated using ( p.lOj ) and produces the last term on the left-hand 
side of ( [^.12| ). The single- length scaling result ~ t^^^ is recovered for the phase ordering 
of the bulk nematic. In Fig. |l| the theoretical results for JF in the the three-dimensional 
0(2) model and the determined in cell-dynamical simulations of the bulk nematic |^ 
are compared. The agreement between the two is excellent. At short-scaled distances has 
the form ( |2.15 ) which is also seen in the simulations and is an indication that string defects 
are the dominant disordering agent in the bulk nematic. 

The theory is now structured so that many well-established phase-ordering results for 
the 0(2) model P,p]^ can be applied to the bulk nematic. In particular, the string line 
density rj is related to the auxiliary field s, whose zeros locate the positions of the strings, 
through 



(3.23) 



where the tangent to the string. 



LU 



Vsi X VS2, 



(3.24) 



points in the direction of positive winding number. The calculation performed in Appendix 
A shows that the average line density of string obeys (rj) ~ Lj'^ ~ for late-times, in 
accord with experiments 0. 



11 



IV. MONOPOLES IN THE BULK NEMATIC 



To address experiments that are designed to produce copious amounts of monopoles at 
the quench |]T^ , a theory for the ordering kinetics of bulk nematics is considered in which 
monopoles appear. The model consists of mapping the director n near a monopole directly 
onto a three-component Gaussian auxiliary field m via n = m. Thus the order-parameter is 



Qap = A{m)[rharhf3 - -6af3] 



(4.1) 



Since the isotropic monopole core can be connected to the nematic minimum along the B = 
line on the potential surface (Fig. a biaxial piece does not appear in (^?T|). Equation 
(O) solves the Euler-Lagrange equation (|2.25|) written in terms of m, 



if the amplitude A satisfies 



V A- 



afS 



:A 



Pa 



3dV 
2dA 



(AO). 



(4.2) 



(4.3) 



A similar result was obtained in |30] for equilibrium. For small m, ( |4.3|) indicates that 



A ~ m while for large m the amplitude A algebraically approaches its ordered value of 
1. The dependence at small m indicates that ( [4.1D describes charge 1 monopoles in the 
nematic |T8|. This is also evident geometrically since fh (and thus n) is a radial vector field 
near the monopole. 

At late-times, using ( [4.1D with A ^ 1 the order-parameter correlation function ( |2.26D is 



(4.4) 



In contrast to the string models considered earlier, the expression (|4.4| ) for the order- 
parameter correlation function in the monopole model is new. The Gaussian average in 



( [4.4| ) is computed in Appendix B. In the late-time scaling regime C(r, t) can be written in 
terms of the scaled length x = r/Lm{t), where Lmit) is the typical monopole separation. 
Thus C{r,t) = T{x) with 



^=l + :^(sin-V-7/) 



(4.5) 



and 7 = 1/a/1 — The auxiliary field correlation function / is defined in ( |2.9|) . The 
scaling function JF satisfies the scaling equation ( |2.12| ) with Lm{t) ~ t^/^. The development 
of this result closely parallels that of the string case considered earlier. The only difference 
between the scaling results for this model and those for the 0(3) model is that the relation 
between JF and / is ( [4.5| ) instead of ( p.ll| ). 

Since m is smooth, / has power series expansion that is analytic at small x. By using 
this expansion in ( |2.12| ) and ( ^.5| ) the small-x behaviour of T is found to be 



:f{x) 



37r 
2/i(i 



x^ + 




(4.6) 
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The non-analytic term in JF, also found in the 0(3) model ( p.l6| ), is due to the presence of 
point monopole defects. Using a fourth-order Runge-Kutta scheme, the non-linear eigenvalue 
problem represented by ( p.l2| ) and is solved for d = 3. The eigenvalue is /i = 1.27306 . . ., 



which differs from the value /i = 0.5558 . . . for the 0(3) model [@]. The function is plotted 
in Fig. ^ along with the scaling function for order-parameter correlations in the three- 
dimensional 0(3) model. Fig. ^ also compares the cell-dynamical simulation data for the 
bulk nematic to the function J-', equation ( [4.51 ). The function JF does not describe the 
simulation data as well as the string model, showing deviations at short distances. These 
deviations are expected since the structure of the theory at short distances ( [4. 61 ) represents 
the wrong defects (monopoles) instead of the correct ones (strings). 

Since the zeros of m locate the monopole cores, the monopole number density n can be 
expressed in terms of the auxiliary field m 



as 



n = 5(m) \det{dimj)\. (4.7) 

where the quantity between the absolute value signs is the Jacobian for the transformation 
from real space coordinates to auxiliary field variables. From the development in the 



average monopole number density obeys (n) ~ ~ This result holds only for 

monopole annihilation in the absence of strings, the case considered in this section. In the 
experiments of Chuang et al. where monopole annihilation occured in the presence of 
strings the average monopole density was observed to decay faster, with (n) ~ t~^. 



V. DISCUSSION 

The dominant scaling behaviour observed during ordering in the bulk nematic is well- 
described by the the theory presented here in which string defects are the major disordering 
agents. The growth law Lg ~ t^/^ is recovered, leading to an average string line density {if) 
that decays as (r/) ~ t~^, as seen in experiments [0. The theoretically determined scaling 
form for order-parameter correlations in the bulk nematic is shown analytically to be exactly 
that for the three-dimensional 0(2) model 0, and this is in excellent agreement with the 
simulation results (Fig. 0). This paper addresses the issue of biaxiality near the string 
cores and demonstrates that it is irrelevant to the leading order scaling properties of the 
system. However, the theory is capable of being extended into the pre-scaling regime, where 
biaxiality may play a role in the dynamics. 

The major accomplishment of this work is the explicit demonstration of the isomorphism 
between the late-stage ordering in the bulk nematic and the late-stage ordering in the three- 
dimensional 0(2) model, within the Gaussian approximation. It is shown that, in the 
scaling regime, the order-parameter equations of motion for the 0(2) model ( p.l|) and the 
bulk nematic (|2.19| ) produce the same scaling equation ( |2.12| ) for the correlation function. 
The essential element in the present theory, which was missing in earlier theories |10|, is 
the mapping ( p.l3| ), which explicitly includes string defects and makes a direct connection 
with the 0(2) model. As a consequence, results for the 0(2) model, such as string and 
vortex density correlations or conservation laws involving string densities 



can 



be directly applied to the bulk nematic. 
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This paper also presents a model for bulk nematics in which monopoles appear. The 
model is applicable to situations where monopole-antimonopole annihilations occur in iso- 
lation from string defects. Such scenarios have been realized experimentally [|I^, and the 
data is suggestive of the growth law ~ t^^"^ predicted by the theory. However, to prop- 
erly treat monopole dynamics in the presence of strings, theories that include interactions 
between string and monopole defects are required. This interesting aspect of the problem is 
under current investigation [ITB]. 
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APPENDIX A: 

This appendix presents the calculation of the average line density for strings, (r^), for the 
0{n) model m d = n + 1 spatial dimensions, defined as 

(ji) = m\^\) (Ai) 

with 



The form (|A2|) , where e is the fully antisymmetric tensor, generalizes the definition (|3.24|) 
for the tangent vector to a string in ci = 3. The one-point average (|A1|) can be written in 
an integral form 



/n+l,n 
n < 1^(01^(0 (A3) 

in terms of 



1 

n\ 

and the one-point reduced probability distribution G{^), given by 



^a(0 = • • • C (A4) 



n+l,n 

G{i) = m n '^(c-v,..)). (A5) 



The Gaussian average in (^3) is straightforward to evaluate by first writing the 5-functions 



in the integral representation and then performing the resulting standard Gaussian integrals. 
One finds 
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n+l,n ( C^\2 



with the definitions 



(27r5o(t))"/2 (27r5(2))"("+i)/' 



(A6) 



5o(t) = -([.-'(0,t)f) 



n+l,n 



n n 



"I" U = \U = \ 



(AT) 
(A8) 



In this theory S^'^^ = l/{n + 1) []T6|. Substitution of ( |A^ ) in ( [A3| ) produces the final form 
for the average fine density of string: 



- gi2) -|V2 

with the n-dependent constant C„ defined as 

1 n+l,n n+l,n 

= ^^^H^mV^ / n <l-(0|exp- E &■ 



(A9) 



(AlO) 



For n = 2 it can be shown that C2 = 1 |[T6|. Since So{t) ~ t at late-times, the average line 
density of string scales like 



-n/2 



(All) 



In particular, for n = 2, {rj) t ^. 



APPENDIX B: 



This appendix outlines the evaluation of the average 

A = ([m(f,t) ■m(0,t)]^) 



(Bl) 



appearing in the correlation function ( [4 .41 ) for the monopole model. For an n-component 
Gaussian m field, the average A can be written in the integral form 



in terms of the two-point reduced probability distribution M 



(B2) 



$(fl,X2) 



7 



271 



7 

exp — — {xl + xl — 2fxi ■ X2) 



(B3) 



where the auxiliary field correlation function / is defined in ( p.9| ) and 7 = 1/ \/l — f^. The 
identity 
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1 1"^ 

= 2 driTi exp — x^r^ 
'1 r -'0 



(xi) 



allows A to be written Gaussian integral 



(B4) 



A = lim / driTi / (ir2r2 



(B5) 



with 



^x{xi,X2,ri,r2) 



Pi' 



SL 

27r 



exp -(r^ + y - + y + 7^-^ ■ X2. (B6) 



The integrals over xi and X2 in ( 1B5| ) are readily done. After differentiating with respect to A 
and setting A = 1, the integral over ri is performed. After a change of variables, y = (^2)^, 
the following integrals remain: 

A = 2-W2+^) r dy{y + ^r\y + V"/^ + 2-^-/'^'Wf H dy{y + trY'iy + \)-^-/'^'\ 
Jo Z Z Jo z z 

(B7) 



These integrals can be expressed in terms of hypergeometric functions F [^], giving 



In ■y 

A = l- — ^ for n = 2 



A = l + 



72 - 1 

in -I) 
72(n-2) 



(F[l,l;n/2;f]-l)-l 



for n > 2. 



(B8) 
(B9) 



In particular, for n = 3, equation ( p9|) gives 



A = l + — 

T 



(7sin"V-/)-l 



(BIO) 



which leads to ( [4.5|) for JF in the nematic with monopoles. 
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FIGURES 




FIG. 1. The scaling function J-{x) for order-parameter correlations as a function of the scaled 
length X = r/Ls{t) for the three dimensional 0(2) model |Q] is represented by the solid line. 
As shown in Section III.B this function exactly describes order-parameter correlations in bulk 
nematics. The cell-dynamical simulation data of Blundell and Bray ||5| for this quantity in a bulk 
nematic are also shown, as circles. The abscissa of the simulation data is scaled so as to give the 
best fit to the theory. 
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FIG. 2. The scaling function J-{x) as a function of the scaled length x = r/Lm{t) for or- 
der-parameter correlations in the theory for monopoles in bulk nematics, discussed in Section IV 
is represented by the lower curve. The upper curve is the scaling function for order-parameter 
correlations in the three-dimensional 0(3) model, from Q. The circles are the cell-dynamical sim- 
ulation data of Blundell and Bray for JF in a bulk nematic. The abscissa of the simulation data 
is scaled so as to give the best fit to the theory for monopoles in bulk nematics. 
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0.0 O.S 1.0 1.5 

A 

FIG. 3. Contour plot of the potential surface V{A, B) (|2.22| ) for the bulk nematic in the physical 
region A > 0, i? > 0. The uniaxial amplitude is A and B is the biaxial amplitude. Darker 
shades indicate lower regions on the potential surface. The isotropic maximum at {A, B) = (0, 0) 
corresponds to a monopole core, the saddle at (^4, B) = (1/4, 3/4) to a string core, and the minimum 
at (j4, = (1,0) to the bulk uniaxial state. 
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